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ABSTRACT 

Recent extreme ultraviolet spectroscopic observations indicate that slow magnetosonic waves are present in active 
region (AR) loops. Some of the spectral data were also interpreted as evidence of fast (~ 100-300 km s _1 ) quasi- 
periodic flows. We have performed three-dimensional magnetohydrodynamic (3D MHD) modeling of a bipolar 
AR that contains impulsively generated waves and flows in coronal loops. The model AR is initiated with a dipole 
magnetic field and gravitationally stratified density, with an upflow-driven steadily or periodically in localized 
regions at the footpoints of magnetic loops. The resulting flows along the magnetic field lines of the AR produce 
higher density loops compared to the surrounding plasma by injection of material into the flux tubes and the 
establishment of siphon flow. We find that the impulsive onset of flows with subsonic speeds result in the excitation 
of damped slow magnetosonic waves that propagate along the loops and coupled nonlinearly driven fast-mode 
waves. The phase speed of the slow magnetosonic waves is close to the coronal sound speed. When the amplitude 
of the driving pulses is increased we find that slow shock- like wave trains are produced. When the upflows are 
driven periodically, undamped oscillations are produced with periods determined by the periodicity of the upflows. 
Based on the results of the 3D MHD model we suggest that the observed slow magnetosonic waves and persistent 
upflows may be produced by the same impulsive events at the bases of ARs. 

Key words: magnetohydrodynamics (MHD) - Sun: activity - Sun: corona - Sun: flares - Sun: oscillations - Sun: 
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1. INTRODUCTION 

Observations of slow-mode oscillations in hot loops 
associated with flows were seen in X-ray emission in the 
past (e.g., Svestka 1994). Observations with the Solar and 
Heliospheric Observatory (SOHO) Solar Ultraviolet Measure- 
ments of Emitted Radiation (SUMER) instrument (e.g., Kliem 
et al. 2002; Wang et al. 2002, 2003a, 2003b) show that stand- 
ing slow magnetosonic waves are excited in hot (~6-9 MK) 
coronal loops by impulsive flare-like phenomena. Recent Hin- 
ode satellite Solar Optical Telescope (SOT) observations, and 
Transition Region and Coronal Explorer (TRACE) satellite ex- 
treme ultraviolet (EUV) imaging observations show that some 
coronal loop oscillations are associated with flows and jets in 
cool (~1 MK) coronal loops (Aschwanden et al. 2002; Ofman 
& Wang 2008). The difficulties in distinguishing between slow 
magnetosonic waves and quasi-periodic flows in some EUV 
time-sequence images of coronal loops are well known (e.g., 
Berghmans & Clette 1999). The use of spectral data in addi- 
tion to stereoscopic EUV imaging to help resolve whether the 
observed event is slow magnetosonic wave or quasi-periodic 
flow is a promising approach (e.g., Marsh et al. 2009; Marsh 
& Walsh 2009). However, the quasi-periodic signatures prop- 
agating along the loops seen in EUV difference images and 
quasi-periodic Doppler blueshift in spectral emission lines are 
still being debated regarding their wave/ flow nature due to the 
limited quality of the present data: e.g., insufficient spectral and 
spatial-temporal resolutions and small signal-to-noise ratios (De 
Pontieu & McIntosh 2010; Verwichte et al. 2010; Wang et al. 
2009, 2010; Tian et al. 2011). 
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Linear MHD models of slow magnetosonic waves in coro- 
nal loops with steady background flows developed in the past 
were limited to straight cylindrical structures and linear waves 
(e.g., Nakariakov & Roberts 1995). Impulsive excitation of 
slow magnetosonic waves without flows in bipolar magnetic 
structures was recently considered in two-dimensional mag- 
netohydrodynamic (MHD; Selwa et al. 2005, 2007), and re- 
cently in three-dimensional (3D) MHD (Ofman & Selwa 2009; 
Selwa & Ofman 2009) models. In these studies it was found 
that the simultaneously excited fast magnetosonic waves could 
facilitate the establishment of standing slow magnetosonic 
waves in curved magnetic field geometry. Recently, Schmidt & 
Ofman (2011) studied slow magnetosonic waves in an expand- 
ing coronal loop associated with a coronal mass ejection us- 
ing the 3D MHD model. The propagation and dissipations of 
slow magnetosonic waves in coronal plumes with background 
solar wind flow were studied using the 2.5D MHD model 
(Ofman et al. 2000). Oscillations of multiple cylindrical coronal 
loop strands that include coupled slow waves, transverse waves, 
and the nonlinear interaction with background steady flow were 
studied using the 3D MHD model (Ofman 2009). 

The present study expands previous models by including for 
the first time impulsive and non-steady flow injection that excites 
the waves in magnetic loops of a bipolar active region (AR) 
in the full 3D MHD model. The nonlinear interactions of the 
oscillations and flows are studied in realistic AR geometry in hot 
coronal AR plasma model stratified by solar gravity. The results 
demonstrate that injection of steady and non-steady flows at 
the footpoints of coronal loops excite compressional waves in 
longitudinal and transverse directions to the loop’s magnetic 
field. The model shows that the same impulsive process at the 
corona-transition region boundary can produce flows and waves 
in coronal AR loops. 
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2. OBSERVATIONAL MOTIVATION 

It was suggested that the standing slow-mode waves in hot 
coronal loops are excited by small (or micro-) flares at one 
footpoint (Wang et al. 2003a, 2005; Wang 2011). Wang et al. 
(2005) examined the evolution of Fe xix and Fe xxi line profiles 
in the initial phase for 54 oscillations and found that nearly half 
show the presence of two spectral components: the stationary 
component and the shifted component that reaches maximal 
Doppler velocity on the order of 100-300 km s _1 . These features 
suggest that the initiation of waves is closely associated with hot 
plasma upflows produced in flares. 

The oscillation event, occurring at 11:20 UT on 2002 April 
17 above the northwest limb, is such a clear example (see 
Figures 1-4 in Wang et al. 2005). This event was first studied by 
Wang et al. (2005) using SOHO / SUMER observations and was 
further analyzed by Wang (2011) using RHESSI data to explore 
the trigger mechanism. SUMER observed two brightenings in 
the Fexix line located at the intersections of the slit with a 
coronal loop. The evolution of Fexix spectral profiles shows 
the emergence of a second component, which has a Doppler 
velocity of up to 230 km s -1 , measured by double Gaussian fits. 
The slow-mode waves were seen to be excited following this 
initial hot flow injection. RHESSI observed a hard X-ray source 
near the northern footpoint of the oscillating loop, suggesting 
that a hot flow pulse produced by a microflare (e.g., via the 
explosive chromospheric evaporation) moves upward along the 
loop, leading to the two brightenings as observed by SUMER. 
This example indicates that the slow-mode wave may be driven 
by the hot flow that is produced by energy release near the loop’s 
footpoint. 

Similarly, the association of continuous upflows and prop- 
agating slow waves observed in coronal loops near the edge 
of ARs by the Hinode EUV Imaging Spectrometer (EIS) sug- 
gests that they may result from the same mechanism, as- 
sociated with tiny nanoflare-like small-scale energy release 
(Wang et al. 2010; Nishizuka & Hara 2011). The SUMER 
oscillation events have high recurrence rate (two to three 
times within 2 hr), consistent with the hypothesis of their 
microflaring origin. We speculate that the nanoflares could 
produce nearly continuous upflows near the loop’s footpoints 
because of their much higher recurrence rate, while larger, 
lower- frequency flares are responsible for the discernible quasi- 
periodic wave disturbances that propagate upward more than 
70 Mm above the footpoints. A noticeable difference is de- 
tected between the upflow-related phenomena observed by EIS 
and SUMER. The EIS continuous upflows and propagating 
slow waves are most clearly observed in warm (1-2 MK) coro- 
nal loops, while those detected by SUMER are observed in 
hot (>6 MK) flaring loops. This observational difference can 
be understood by recent nanoflare models (see the review by 
Klimchuk 2009) that predict the presence of a small amount 
of very hot (>5 MK) plasma. There are two main reasons that 
make hot spectral lines hard to detect: the shorter lifetime and 
the lower densities of the hot material compared to the warm 
plasma. Therefore, we model these phenomena in hot coronal 
loops as a first reasonable step with the aim of understanding the 
relationship between flows and waves excitations by impulsive 
events at the loops’ footpoints. 

3. NUMERICAL MODEL 

Here, we describe the numerical 3D MHD model used to 

study the hot AR loops with flows and waves. The resistive 
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Figure 1. Initial dipole magnetic field (white curves) used for the model AR. 
The intensity scale shows the magnetic field magnitude at the base of the AR 
and the white arrows on some field lines show the direction of the magnetic 
field. The dimensions are in normalized units. The arrow shows the location and 
direction of the injected flow. 

(Animations and a color version of this figure are available in the online journal.) 

3D MHD equations are solved with the gravity and isothermal 
energy equation on a Cartesian 25 8 3 grid using the modified 
Lax-Wendroff method with fourth- order stabilization (smooth- 
ing) term. The initial magnetic field is a dipole (see Figure 1) 
with gravitationally stratified normalized density given by 

p = p^n/ao+z-z^-o.i ]/h (1) 

where H = 2kBToR s /(10GM s m p ) is the normalized gravita- 
tional scale height, R s is the solar radius, k B is Boltzmann’s 
constant, To is the background (isothermal) temperature, G 
is the universal gravitational constant, M s is the solar mass, 
m p is the proton mass, and p 0 = 1. The following nor- 
malization parameters were used in the present study: dis- 
tances in units of a = R s / 10, magnetic filed Bo = 100 G, 
T 0 = 6.3 MK (to model the hot coronal AR loops), number 
density no = 1.38 x 10 9 cm -3 , resulting in Alfven speed Vao = 
Bo/(4irno) l/2 = 5872 km s -1 used to normalize the velocities, 
the Alfven time z A — a/ V A o = 11.9s used for time normaliza- 
tion, the isothermal sound speed C s = 323 km s -1 = 0.055 V A o, 
and the corresponding plasma fio = (2C 2 )/(V^ 0 ) = 0.00605. 
The model is based on the 3D MHD model of a bipolar AR de- 
veloped initially by Ofman & Thompson (2002). This model and 
its extensions and variants were used successfully in a number of 
studies of waves in ARs and in coronal loops (e.g., Ofman 2007; 
McLaughlin & Ofman 2008; Selwa & Ofman 2009, 2010; Selwa 
et al. 201 la, 201 lb). As in previous studies we set the Lundquist 
number to S = 10 4 * , i.e., the resistive diffusion timescale is four 
orders of magnitude longer than the Alfven time. Since we are 
dealing with timescales O(100z A ) — the resistive dissipation is 
insignificant. Numerical convergence study has shown that the 
results reported below are not affected by numerical dissipation. 

In Figure 2 the initial density, p, the fast magnetosonic 
speed for perpendicular propagation Vf = (V A + C 2 ) 1/2 , where 
V A = B(x, y, z)/p(x, y, z) x ^ 2 V A o is the local Alfven speed, 
and the plasma /3 in the a z -plane at y = 0 of the model 
AR are shown. The gravitational stratification of the density 
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Figure 2. Initial density (left), fast magnetosonic speed (middle), and plasma /3 (right) in the xz-plane at t = 0 in the model AR. The contours on Vf show the 50%, 
25%, and 12.5% levels of the maximal value, and the contour /3 = 1 shows the heights where this value is reached. 

(A color version of this figure is available in the online journal.) 


is evident, and the increase of Vf in the lower central region 
due to the bipolar magnetic field of the AR is apparent. The 
plasma ft increases with height due to the rapid decrease of 
the dipole magnetic pressure ~B 2 compared to the decrease 
of the thermal pressure with height. 

Non-steady and impulsive flow injection along the field is 
introduced at the lower boundary of the model AR as 

V = V 0 (x, y,z = Zmin, t)B/\B\, (2) 

where 

Vo(x, y,Z = Zmin, 0 



with 

A v (t) = to < t < t max , (4) 

for impulsively excited continuous flow or 

A v (t) = A 0 (l — cos cot)/2, t > to (5) 

for periodic upflow. Here, the parameters are z m in = 1, xo = 0.8 
for the short magnetic loop (length = 2.13), xo = 1.2 for the 
long magnetic loop (length = 4.53), yo = 0, wo = 0.12, 
co = 0.21, to = 0, and the values A 0 = 0.01 and 0.05 
that correspond to subsonic upflow were used. The additional 
power of two on the right-hand side of Equation (3) models 
an upflow with sharper- than- Gaussian cross-sectional profile. 
The upflow is imposed at the lower boundary z = Zmin = 1 
at a circular area with diameter 2wo centered at the right 
footpoint of the loops. Downflows through the lower boundary 
are allowed at the other footpoint of the loops. The rest 
of the boundary conditions at z = 1 are v x = v y = 0 
outside the up flow/downflow regions, with fixed B, density p, 
and v z extrapolated from the interior points. These boundary 
conditions approximate the line-tied boundary conditions for 
coronal loops (e.g., Lionello et al. 1998; McLaughlin & Ofman 
2008), modified allowing inflow and outflow. The line-tied 
boundary conditions approximate the effects of the dense 
photosphere on the loop’s magnetic footpoints, and ensure the 
wave reflection necessary for establishment of standing waves 


with nodes at the footpoints. Open boundary conditions are used 
at the top and side planes of the computational box, allowing 
wave propagation outward through the boundary with negligible 
reflection. 

The magnetic field lines of the loops and the dependence of 
the loops parameters on the coordinate along each loop, s, are 
shown in Figure 3. It is interesting to note that the normalized 
magnetic field strength B and the fast magnetosnic speed Vf 
along the loop decrease with height by a factor of three for 
the short loop and by an order of magnitude for the long 
loop. The significant variation of the magnetic field and the 
fast magnetosonic speed is usually neglected in linear coronal 
seismology models (see the review by Nakariakov & Verwichte 
2005). The density along the loops varies by ~10% for the 
short loop and ~23% for the long loop. The small variation of 
the density with height (compared to cool 1 MK loops) is the 
result of the high temperature (6.3 MK) and the low height of 
these loops. 

4. NUMERICAL RESULTS 

In Figures 4-14, the results of the 3D MHD model of steady 
and periodic inflows and associated waves are shown for the 
initial state and parameters described in the previous section. 
The evolution of the flows and oscillation is also evident in 
the animations accompanying the online version of the journal. 
First, we report the results for the steady inflow cases. In 
Figure 4, the snapshots of the density and the velocity in the 
xz-plane at the end of the run ( t = 280 ta) of steady inflow 
with Vo = 0-05 are shown. The intensity scale indicates the 
magnitudes of the density and the velocity, and the arrows 
in the velocity images indicate the direction of the flow for 
velocities greater than 6% of the maximal value. The short loop 
results are shown in Figures 4(a) and (b), and the long loop 
results are in Figures 4(c) and (d). It is evident that the steady 
inflow leads to the formation of a steady short loop in the low ~/3 
region of the model AR. The short loop is filled with inflowing 
coronal plasma and reaches higher (nonuniform) density than 
the surrounding plasma by a factor ~1.5. Since the EUV line 
emission in optically thin coronal plasma is proportional to 
n 2 this loop would appear twice brighter than the same loop 
preceding the inflow in background- subtracted EUV images. In 
Figure 6, the cut in the xy-plane at z = 1.26 of the density 
and the velocity of the AR with the short loop are shown. The 
locations of the upflow and downflow that correspond to the 
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Figure 3. Initial structure of the long (solid line style) and short (dashed) loops, (a) The magnetic field lines at the center of the loops in the xz-plane at y = 0. The 
red circles A, B, C, D mark the locations of temporal evolution plots of the variables shown below, (b) The magnetic field intensity along the loops, where 5 is the 
coordinate along each loop, (c) The density along the loops, (d) The fast magnetosonic speed, V f(s) along the loops. 

(A color version of this figure is available in the online journal.) 
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Figure 4. Snapshots of the density and the velocity in the xz-plane at t = 280ta with steady inflow (Vo = 0.05). The intensity scale shows the magnitude of the density 
and the velocity, and the arrows indicate the direction of the flow for velocities greater than 6% of the maximal value, (a) Density for the short loop, (b) velocity in the 
short loop, (c) density for the long loop, (d) velocity in the long loop. The rectangles marked E-H mark the locations of the time-distance plots shown in Figure 5. 

(A color version of this figure is available in the online journal.) 
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Figure 5. Time-distance plots for the boxed regions E-G of the short loop and H of the long loop shown in Figure 4. Panel E is the time-distance plot obtained from 
the cut E in the z-direction of the short loop centered at x = 0. Panel F is the time-distance plot of the cut F in the x-direction centered at z = 1.82. Panel G is the 
time-distance plot of the cut G in the x-direction centered at z = 1.27 of the short loop. Panel H is the time-distance plot of the cut H in the x-direction centered at 
z = 1.82. In panel H, the initial transient time is removed to show better the subsequent evolution of the loop’s density. In all panels the temporal cadence is 5ta- Both 
displacement and density oscillations are present showing evidence for the transverse and compressional waves in the x-z plane of the loops. 

(A color version of this figure is available in the online journal.) 


loop with higher than the surrounding density are evident. In 
addition the circular structures in the density and the velocity 
are formed due to the compression and rarefaction produced by 
the flow in the AR and the fast magnetosonic waves emitted in 
the initial impulsive stage of the loop formation. 

The time-distance plots shown in Figure 5 for the boxed 
regions E-H marked in Figure 4 demonstrate the details of 
the loop oscillations. The widths of the boxes are seven gird 
points in space and the time-distance plots are smoothed with 
15 points running average. The snapshots of the boxes are 
taken every 5 Ta. The time-distance plot E is aligned along the 
z -direction centered at x = 0 through the apex of the short 
loop. It is evident that the temporal oscillation of the loop in the 
z-direction is small, and the density variations are due primarily 
to compressive density oscillation, following an initial transient 
stage (the time for the inflow to reach the apex of the loop). The 
time-distance plot F is aligned along the x-direction centered 
at height z = 1.55 through the apex of the loop and shows 
the oscillatory displacement of the loop top, with evidence of 
compressive density oscillations. By comparing time-distance 
plots E and F it is evident that the dominant transverse (kink) 
oscillation of the loop is in the x-direction (see, Figure 8 below). 


The time-distance plots G (centered at z = 1.27) and H 
(centered at z = 1.82) are cut in the x-direction through the 
footpoints of the short and long loops, respectively. Note that 
the positive values of x correspond to the right side of the plots 
in Figure 4. Here, the transverse oscillations of the footpoints, 
as well as the effects of density compression are evident. The 
displacements of the loops’ two footpoints are in-phase, i.e., 
the loop sways left and right in the loops’ curvature plane, 
whereas the intensity oscillations near the two footpoints are 
in anti-phase, a signature of the fundamental slow-mode wave, 
which has the density anti-nodes at the footpoints. These loop 
oscillations are similar to the second harmonic of the vertical 
kink mode discussed by Wang et al. (2008) in connection with 
TRACE EUV coronal AR loop oscillations. 

Figure 7 shows the temporal evolution of the velocity com- 
ponents, the perturbed magnetic field components, and the per- 
turbed density at the right footpoint of the short loop (point 
A = (0.59, —0.01, 1.27)) with steady inflow. At this location 
the velocity along the loop has two components in the Carte- 
sian frame (x and z) as evident in Figure 3(a). It is evident 
that it takes about 15 Za for the perturbation along the loop to 
reach this position. The loop contains both compressive and 
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Figure 6. Same as Figure 4, but in the xy - plane at z = 1.26 for short loop, (a) Density, (b) velocity. 

(An animation and a color version of this figure are available in the online journal.) 



Figure 7. Temporal evolution in a short loop of the velocity components (V*: dots; V y : dashes; V z : solid) (top panel), the perturbed magnetic field components (B x : 
dots; B y \ dashes; B z : solid; | A B | : dot-dashes) (middle panel), and the perturbed density (lower panel) at the right footpoint position A = (0.59, — 0.01 , 1 .27) for steady 
inflow. The red dashed lines through the third peak of V z and the peak of | A B | help evaluate the phase shifts between these variables. 

(A color version of this figure is available in the online journal.) 
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Figure 8. Temporal evolution in a short loop of the velocity components ( V x : dots; V y : dashes; V z : solid) (top panel), the perturbed magnetic field components (B x : 
dots; B y \ dashes; B z : solid; |AZ?|: dot-dashes) (middle panel), and the perturbed density (lower panel) at the apex of the loop B = (0.01, —0.01, 1.6) for steady inflow. 
The red dashed line marks three peaks of An helping evaluate the phase shift between the variables. 

(A color version of this figure is available in the online journal.) 


transverse oscillations, as evident in the time dependence of the 
velocity and magnetic field component oscillations, as well as 
the density fluctuations. The transverse waves are excited by 
the momentum of the injected velocity pulse and the centrifugal 
force exerted by the inflow on the curved (dipole) field lines, 
with the Lorentz force acting as the restoring force. The slow- 
mode wave period is 65 Ta (determined from the last oscillation 
period of An), and the oscillations damp in several periods due 
to leakage at the footpoints and outside the loop as a result 
of finite ft allowing oblique propagation and leakage of the 
slow magnetosonic waves in curved magnetic geometry. The 
quarter period phase shift between the density and the velocity 
components along the loop indicate the presence of standing 
slow-mode waves is established quickly. The presence of the 
slow-mode wave is also evident from the nearly anti-phase re- 
lation between |A2?| = (A + A + A and An (see the 
vertical red dashes lines in Figure 7). The anti-phase relation 
is not exact due to the effects of nonlinearity and non-modality 
(i.e., the oscillation are not exact normal modes of the loop) of 
the oscillations. The effects of nonlinearity are evident in the 
non-sinusoidality of the oscillations at the initial stages of the 
evolution. The calculations were repeated with smaller velocity 
amplitude (Vo = 0.01) confirming that these effects are due to 
nonlinearity and the magnitude of the initial pulse. 


Figure 8 shows the temporal evolution of the velocity com- 
ponents, the perturbed magnetic field components, and the 
perturbed density at the apex of the short loop (point B = 
(0.01, —0.01, 1.60)) with steady inflow. The red dashed line 
marks three consecutive peaks of the density perturbation, An. 
It is evident that A B x is in-phase with An indicating the presence 
of fast-mode wave and the corresponding kink oscillation of the 
curved loop in the vz- plane induced at the impulsive initial stage 
of the inflow. The variations in An and A B x are caused by the os- 
cillations of the nonuniform loop with respect to the fixed point 
B , with contribution from compressive effects (see Figure 5). 
The period of the oscillation is 35 Ta based on the second and 
third peaks. This period is shorter than the slow wave period and 
longer than the kink mode period of the fundamental mode for a 
cylinder of length L (e.g., Roberts et al. 1984) given by 2 L/Ck, 
where C* = Va(2/(1 + p 0 /pi) 1 ^ 2 = 1.1 Va, with the ratio of the 
density inside and outside the loop Pi/p 0 = 1.25/0.78 = 1.5 
from the model calculations. Noting that V a decrease with height 
in the loop we get Ck = 0.35 V^o at the apex of the loop. How- 
ever, the theoretical kink mode period of 12 ta is a factor of 
three shorter than the period found in the model. Thus, the ob- 
served fast mode is not a linear fundamental kink mode of the 
loop but nonlinearly driven fast-mode oscillation with similar 
properties to the second harmonic of the vertical kink mode. 
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Figure 9. Temporal evolution in a long loop of the velocity components ( V x : dots; V y : dashes; V z : solid) (top panel), the perturbed magnetic field components ( B x : 
dots; By : dashes; B z : solid; | A B | : dot-dashes) (middle panel), and the perturbed density (lower panel) at the right footpoint position C = (1.08, —0.01, 1.82) for steady 
inflow. The red dashed lines through the third peak of V z and the peak of | A B | help evaluate the phase shifts between these variables. 

(A color version of this figure is available in the online journal.) 


Further justification of this interpretation is obtained by in- 
specting the phase and period relation between A B x , An, and 
AV X . It is evident that AV X period is twice as long as the A B x , 
An, and peaks, nodes, and minima of AV X all correspond to 
peaks of An and A B x , suggesting that the fast mode is driven 
nonlinearly by a term that depends on (AV X ) 2 . Such a second- 
order term appears in the v component of the momentum equa- 
tion ~pV x (dV x )/(dx ) = (p/ 2)(d(V x ) 2 )/(dx). The nonlinearly 
driven fast-mode waves damp as the driving amplitude of the 
slow waves decreases with time. 

Figure 9 shows the temporal evolution of the velocity com- 
ponents, magnetic field components, and the perturbed den- 
sity at the right footpoint of the long loop (point C = 
(1.08, —0.01, 1.8)) with steady inflow. The long loop exhibits 
similar evolution and contains both longitudinal and fast-mode 
waves induced by the flow. The slow wave period is ~ 136 ta, 
determined from the last half period of An oscillation. The pres- 
ence of the slow-mode wave is also evident from nearly anti- 
phase relation between | A B | and An (see the vertical red dashes 
lines in Figure 9). The anti-phase relation is not exact due to the 
effects of nonlinearity and non-modality (i.e., the oscillation are 
not exact normal modes of the loop) of the oscillations. As in 
the short loop case the waves damp due to leakage of the slow 


wave outside the loop. Evidence for the formation of the slow 
standing wave is seen in the nearly quarter period phase shift 
between the velocity along the loop and An. However, in the 
case of the long loop the periods of the oscillations is longer 
as expected, and the formation of the slow standing wave takes 
longer time than in the short loop. By comparing the short and 
long loop oscillations it is evident that the periods of the os- 
cillations and the wavelengths are related to the loop lengths. 
The ratio of their periods is about a factor of two, close to the 
loop’s length ratio of 2.15. Wave period values in both loops are 
in qualitative agreement with the theoretical expression 2 L / C s 
for the fundamental mode of the slow magnetosonic wave in a 
straight cylinder (e.g., Roberts et al. 1984) with the short loop 
in better agreement. However, quantitatively the slow wave pe- 
riods are smaller by 16% for the short loop and by 17% for 
the long loop case. The difference is mainly due to the curved 
magnetic geometry and nonuniformity of the loops. 

In the following, we report the results for the periodic flow 
injection cases. The density and velocity distributions in the 
xz - plane (at y = —0.0136) cuts of the AR produced by periodic 
flow injection at the lower boundary at (0.8, 0, 1) at times 
t = 22.5, 105, 188 Ta for the short loop are shown in Figure 10. 
The loop density structure shows time-dependent variation due 
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Figure 10. Density (top panels) and the velocity (lower panels) in the xz-plane (y = —0.0136) at t = 22.5, 105, 188 ta (from left to right) in a short loop injected 
periodically with flow. The arrows show the direction of the flow and the intensity shows the magnitude of the flow and the density. 

(An animation and a color version of this figure are available in the online journal.) 



Figure 11. Temporal evolution of the velocity components ( V x : dots; V y : dashes; V z : solid), the perturbed magnetic field components (B x : dots; B y : dashes; B z \ solid; 
| A B | : dot-dashes), and the perturbed density at the right footpoint position A = (0.59, — 0.01 , 1.27) of the short loop with periodic inflow. The red dashed line through 
the fourth peak of V z line helps evaluate the phase shift between the variables. 

(A color version of this figure is available in the online journal.) 
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Figure 12. Density (top panels) and the velocity (lower panels) in the xz-plane at t = 22.5, 105, 188 xa (from left to right) in a long loop injected periodically with 
the flow. The arrows show the direction of the flow and the intensity shows the magnitude of the flow and the density. 

(An animation and a color version of this figure are available in the online journal.) 


to the combined effects of the periodic flow injection and 
the compression of AR plasma by the velocity pulse. The 
velocity magnitude and the direction vary approximately in- 
phase with the density showing longitudinal variations along the 
expanding magnetic field of the loop. The temporal evolution of 
the velocity components, perturbed magnetic field components, 
and perturbed density are shown in the short loop at the right 
footpoint at point A = (0.59, —0.01, 1.27) for the periodic 
flow injection in Figure 11. It is evident that after the short 
initial phase as the perturbations reach the above position, 
the loop’s magnetic field, velocity, and density exhibits quasi- 
periodic oscillations induced by the periodic flow injection at 
the footpoint. The velocity fluctuations along the loop are in- 
phase with the density perturbation, and in anti-phase with the 
magnetic field components. The temporal evolution of \AB \ 
shows two dominant frequencies due to the combined effect of 
the driver and the slow-mode wave in the loop. No temporal 
damping or growth of the oscillations is evident in the driven 
case at the steady state, since the fluctuation are injected at 
the right footpoint continuously and exit the loop at the left 
footpoint boundary. The non-sinusoidality of the oscillations is 
evident and is due to the nonlinear interaction between the flow 
periodicity, the resulting variation of the local loop Alfven speed 
that affects the properties of the Alfvenic fluctuations, and the 
waves produced in the loop. The nonlinear effects are reduced 
when the same calculation is repeated with smaller injected 
velocity amplitude (Vo = 0.01 , not shown). 

In Figure 12, the cuts in the xz-plane of the AR’s density 
and velocity injected with periodic flow at the lower boundary 
at (1.2, 0, 1) at times t = 22.5, 105, 188 t a for the long loop 
are shown. It is evident that in this case the flow is divergent as 


expected in the bipolar magnetic field geometry. In the case of 
the long magnetic loop, the density structure is more complex 
than in the case of the short loop that is entirely in the low ~/3 
region of the AR as shown above. Since the long loop reaches 
higher altitude than the short loop, nonlinear steepening of the 
density perturbation fronts are produced due to the decrease 
of the gravitationally stratified density and the increase of the 
plasma with height (see Figure 2). Coupled slow and fast 
magnetosonic shock-like waves are produced as evident in 
the compression fronts in the images, and their direction of 
propagation that is oblique to the direction of the magnetic 
field. Here as well no temporal damping of the oscillations is 
evident. Figure 13 shows the temporal evolution near the apex of 
the long magnetic loop (point D = (0.01, —0.01, 2.5)) for the 
periodically driven upflow case. The shock- like structures are 
evident in the temporal evolution of the velocity components, 
the density, and the magnetic field exhibiting strong nonlinear 
steepening. Here as well, the temporal evolution of \AB\ shows 
two dominant frequencies due to the combined effect of the 
driver and the slow-mode wave in the loop. The structures 
are shock- like (more diffuse compared to shocks), since they 
travel less than a single scale height, and their amplitude is 
significant but not large compared to the sound speed and fast 
magnetosonic speed. Convergence test shows that numerical 
dissipation does not affect significantly the shock-like waves in 
the present calculation. The nonlinear steepening of the flow 
fronts in the loop is due to inflow amplitude comparable to the 
sound speed at the footpoint, and the density stratification in 
long loops (see Figure 3). 

Finally, we summarize the results of the four cases discussed 
above using 3D isocontours snapshots of the perturbed density 
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Figure 13. Temporal evolution of the velocity components (V*: dots; V y : dashes; V z : solid), the perturbed magnetic held components (B x : dots; B y : dashes; B z : solid; 
|AZ?|: dot-dashes), and the perturbed density near the apex of the loop at position D = (0.01, —0.01, 2.5) of the long loop with periodic inflow. The red dashed line 
though the third peak of V x helps evaluate the phase shift between the variables. 

(A color version of this figure is available in the online journal.) 


in Figure 14. The four cases discussed above are: (1) steady 
inflow in the short magnetic loop, (2) periodic inflow in the short 
magnetic loop, (3) steady inflow in the long magnetic loop, and 
(4) periodic inflow in the long magnetic loop. By comparing 
between the isocontour snapshots and animations the effects of 
the magnetic loop sizes and the types of inflow become apparent. 
A steady loop is formed following steady inflow and induced 
damped oscillations in the initial transitory stage (longitudinal 
and transverse) in the short loop that lays entirely in the low-/? 
region of the AR. When the periodic flow is injected in such 
magnetic loop the density structure in the loop is dominated by 
the structure of the flow with small effects due to waves in the 
loop. In the long loop the steady state is not reached, and the 
density of the loop is variable throughout the evolution for both 
types of inflow. When the inflow is periodic the effects of the 
waves induced by the flow are significant in the long loop (for 
the magnitude of the flow in the present study) and the density 
structure exhibits the coupled fast magnetosonic and the slow 
wave shock-like structure in the p ~ 1 region of the model AR. 

5. DISCUSSION AND CONCLUSIONS 

High-resolution and high cadence observations indicate that 
some coronal loop oscillations are associated with flows along 
the loops and with jets at the loops’ footpoints. Whether propa- 


gating intensity disturbances seen in EUV imaging observations 
of the coronal loop structures are flows or waves is still under 
intensive debate. On the one hand the propagating disturbances 
have been interpreted as signatures of slow magnetoacoustic 
waves since their propagating speeds are close to the sound 
speed in the corona, and three-minute or five-minute period- 
icities suggest their association with the leakage of the pho- 
tospheric /?-modes (e.g., Nakariakov et al. 2000; De Moortel 
et al. 2002). Combined stereoscopic EUV and spectral observa- 
tions provide further evidence for slow magnetosonic waves in 
coronal loops (e.g., Marsh et al. 2009; Marsh & Walsh 2009). 
Observations of propagating disturbances in coronal plumes 
were interpreted as slow magnetosonic waves not only based 
on their speed but also on their amplitudes’ dependence with 
height, found to be in agreement with the expected variation 
of the wave amplitude with height in gravitationally stratified 
plasma (Ofman et al. 1999). On the other hand, some recent 
studies suggested that propagating EUV disturbances are most 
likely intermittent transient upflows because of their association 
with persistent outflows at the loop base (e.g., Hara et al. 2008), 
and their correlation with blueward asymmetries of line profiles 
using Hinode / EIS observations (McIntosh & De Pontieu 2009; 
Tian et al. 2011). However, these observational signatures were 
found to be consistent with slow-mode waves as well (Verwichte 
et al. 2010). 
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Figure 14. Perturbed density isocontours due to inflow in the short loop: (a) steady inflow, isocontour level at p = 0.12; (b) periodic inflow, isocontour level at 
p = 0.1. Same for the long loop: (c) steady inflow, isocontour level at p = 0.035; (d) periodic inflow, isocontour level at p = 0.05. 

(Animations of this figure are available in the online journal.) 


Our 3D MHD model results of flows and associated waves in 
hot coronal loops show the close connection between the two 
phenomena and suggest that both may be present simultaneously 
and produced by the same impulsive events at the footpoints 
(or corona/transition region interface) of coronal loops. In 
particular, we study the effects of localized inflow at the 
lower coronal boundary in a bipolar magnetic field of hot 
AR coronal loops. The 3D MHD model reproduced the basic 
physical coupling between steady and periodic flows, waves, and 
loop density and magnetic structures in realistic AR geometry, 
with applications to coronal seismology, extending previous 
studies that use more simplified models and help addressing the 
controversy of flows versus waves discussed above. The main 
limitations of our AR model are the use of idealized magnetic 
structure with dipole field and gravitationally stratified density 
in the initial state. In addition the isothermal energy equation is 
used, simplifying the calculations considerably at the expense of 
calculating the effects of heat conduction (the dissipation of slow 
magnetosonic waves in hot coronal loops by heat conduction 
was studied by Ofman & Wang 2002), heating and cooling of 
the material injected into the loops, and the possible heating of 
the loops by the energy released in the impulsive events. 

We investigate the effects of steadily injected inflow in short 
and long magnetic loops and the associated oscillations in the 
initial impulsive stage. We find that the inflow at the footpoints 


of coronal magnetic loops leads to formation of higher density 
than the surroundings loops due to the siphon effect in the 
curved bipolar magnetic structure. The inflows excite coupled 
longitudinal and transverse oscillations in the loops in the initial 
(impulsive) phase that later damp primarily through leakage 
as the loops reach the new steady state with background flow. 
The fundamental standing slow-mode waves are seen to form 
in short and long loops in agreement with past observations 
and models discussed above. The fast magnetosonic waves 
induced by the inflow and nonlinearity are detected as well. 
In short loops with steady inflow in low ~/3 regions of the AR the 
properties of the oscillations (period, wavelength, and damping 
rate), depend on loop length, and curvature of the magnetic 
structure in qualitative agreement with previous studies of waves 
in loops without background flow (McLaughlin & Ofman 2008; 
Selwa & Ofman 2009, 2010; Selwa et al. 2011a, 2011b). We 
also find that the flows produce impulsively nonlinear fast-mode 
waves in the loops. 

We investigate the effects of quasi-periodic flow injection 
in the localized region at lower coronal boundary and find that 
driven undamped oscillations are produced in the magnetic loop 
of the bipolar AR’s. As a result, the average density of the 
loops increases as well in this case. The period, amplitude, and 
damping rate of the driven oscillations are determined primarily 
by the properties of the driving periodic flows such as their 
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period and amplitude, and the nonlinear interaction with the 
loop structure. In hot and long loops the quasi-periodic inflows 
with subsonic magnitude at the lower boundary can produce 
shock-like nonlinear compressive disturbances that propagate 
obliquely to the magnetic loop direction. 

We conclude that a way to determine whether the observed 
loop oscillations result from slow waves or quasi-periodic flows 
is to evaluate the relation between the oscillation properties 
such as period, phase speed, phase relations, damping time, and 
the various loop parameters such as loop length, density, tem- 
perature, and magnetic field structure, since these relations are 
different for waves and for quasi-periodic flows as demonstrated 
in our study. The presence of the temporal damping of the oscil- 
lations on the timescale shorter than the decay time of the flow 
(as in the case of steady inflow) supports the wave interpreta- 
tion of the observed oscillations. In the case of periodic inflow 
studied here the oscillations do not exhibit damping. We con- 
clude that the excitation of slow-mode (and some transverse) 
oscillations and flows observed with SUMER, TRACE , Hinode , 
and recently with Solar Dynamics Observatory in coronal loops 
may result from the impulsive flow injection of plasma at the 
corona-chromosphere interface at the footpoints of the coronal 
AR loops. Thus, our study suggests that flows and waves in 
coronal loops may result of the same impulsive events that drive 
these phenomena. 
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